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ELEMENTS OF ACTIVE VIBRATION CONTROL FOR ROTATING MACHINERY 


Heinz Ulbrlch* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

Demands for higher frequencies, Improved reliability, reduced noise, and 
Increased longevity, along with safety concerns, require an effectively con- 
trolled rotor dynamics system. The desired requirements are often not met by 
using only passive damping elements; further Improvements can be achieved only 
with the aid of active control. Showing how to construct a controller and how 
to apply active control on rotating machinery Is the intent of this report. 

Different types of actuator concepts and their models are discussed. 

The modeling of rotor systems Is given In such a way that Inclusion of active 
elements can easily be achieved. The chosen modeling procedure, which is both 
physically clear and handy for computer-oriented representation, permits both 
the consideration of all mass and gyroscopic effects and a modular construc- 
tion of the system. The meaning of controllability, observability, and spill- 
over, with regard to actively controlled systems, and a method to check on 
the system properties are discussed. A short Introduction about control con- 
cepts and the optimization of the controller Is given. In the last section 
several real applications demonstrate the design and application of active 
vibration control . 


INTRODUCTION 

In the growing field of rotor system dynamics control, many different 
methods have been applied. The methods that have been proposed for rigid and 
flexible rotors differ In control concepts (refs. 1 to 8). The type of con- 
troller that is most effective depends on the objectives. The controller nor- 
mally makes use of four elements: 

(1) Displacement (proportional) feedback (P part), which allows one to 
change the stiffness of the system or. In other words, to shift the 
elgenfrequencles In the desired manner 

(2) Velocity (derivative) feedback (D part), which permits one to increase 
the system damping to (a) Improve the stability of the system, (b) sta- 
bilize unstable modes of the system, or (c) reduce resonance amplitudes 
that can occur when rotors pass through critical speeds 

(3) Acceleration feedback (A part), which can be used to eliminate or 
reduce the Influence of the mass or Inertia (e.g., Improving the con- 
trollability In relation to applying active forces via the bearing hous- 
ings (ref. 4), or compensating for the Influence of fluid Inertia In 
relation to hydraulic actuators (ref. 2)) 
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(4) Integral feedback (I part), which Is necessary for holding a rotor 
position very precisely (no radial displacement) at a given location, 
or for adjusting bearings to minimize rotor stress (minimum stress Is 
achieved by changing the features of the bearing or the foundation of 
the bearing supports). 

The type of feedback control that Is chosen will depend on one's objec- 
tives and the system Itself. 

To realize active control techniques for Improving the dynamic behavior 
of rotor systems, the following are required: (1) suitable actuators and 

knowledge of their transfer characteristics; (2) modeling of the complete sys- 
tem (including the active components); (3) choosing locations for actuators 
and sensors (affects controllability, observability, and spillover effects); 
and (4) having problem-adapted control methods for hardware as well as 
software. 

The key to successful active vibration control in rotor systems (l.e., 
systems with rotating machinery) lies In the availability of suitable actua- 
tors. Some types of actuator that have been applied In Industry or In the 
test stage In the laboratory will be Introduced In this report. 


ACTUATORS, MODELS, AND TRANSFER CHARACTERISTICS 

In general, there are two different ways to apply control forces on 
rotating shafts In turbine engines, power generators, gearboxes, machine 
tools, and such. The forces can act directly on the rotor or they may be 
applied via the bearing housings. For practical applications the following 
characteristics are very Important: 

(1) Actuators should not only be compact but also be capable of generat- 
ing large forces. 

(2) The amplitude range of the actuator should be at least as high as the 
vibration amplitudes that have to be Influenced. 

(3) The frequency range of the actuator (bandwidth) that can be obtained 
will determine the applications. This means the design strategies for the 
actuators will have to focus to a high degree on the actuator dynamics . In 
this section three different types of actuators will be Introduced and their 
mathematical descriptions of the transfer characteristic will be given. 


Magnetic Actuators 

Magnetic actuators can be divided into two types: magnetic bearings, 

which act directly on the rotor without contact (refs. 1, 2, 9, and 10), and 
electromagnetic actuators, which apply forces to the rotor Indirectly via the 
bearings (refs. 3, 5, and 11 to 13). 

Active magnetic bearing . - A typical assembly of electromagnets (EMX1, 
EMX2, EMY1, EMY2) is represented In figure 1. It Is noteworthy that with this 
arrangement of the magnets, radial movement of the rotor In the y-dlrectlon 
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causes no change In the total air gap of the magnetic circuits irt the 
x-dlrectlon (comparative effect). This Is necessary for decoupl 1 ng the forces 
of the bearing between the x- and y-dl rections . The differential arrange- 
ment of the electromagnets leads to a linear input/output characteristic 
between the forces acting directly on the rotor and the control voltages input 
to the power amplifiers. The power amplifiers have to work, In this case, as 
current sources. This relationship between the Input voltages u x and u y 
and the resulting forces F x and F y is given by 



where k u represents the force-voltage factor; k, the force-displacement fac- 
tor (negative of stiffness); and x and y, the radial displacements of the 
rotor in the respective directions. Figure 2 shows a schematic of a 
control led-orblt magnetic bearing test rig to determine k u , k, and the fre- 
quency characteristic. The rotor Is fixed by ball bearings and Its position 
Is measured by built-in noncontacting displacement transducers. Forces on the 
rotor are simultaneously determined In the x- and y-dlrectlons by quartz 
load washers. In order to decouple the forces Into the two measurement direc- 
tions, special axial linear roller guides are used. This means the quartz 
load washers are attached to the foundation in such a way that a force In the 
x-dlrectlon (and similarly, In the y-directlon) does not influence the force 
signal In the y-directlon (x-dlrectlon); see references 10 and 14. 

Flectromaanetic actuator. - An electromagnetic actuator capable of 

applying forces to the rotor Indirectly via the bearing housing Is shown In 

figure 3. The scheme shows a cross section of this active element, consisting 
mainly of two annular electromagnets acting In differential prlnc pie on a 
pull disk. If the power ampl if Vers work as current sources, the input/output 
characteristic is described by the same equation as that for the magnetic bear- 
ing (eq. (2). The transfer characteristic of the actuator Is shown in fig- 
ure 4(a). The control force is displayed as a function of the control 
current. The deviation of the pull disk from Its neutral position functions 
as a parameter. Figure 4(b) gives an Impression of the frequency characteris- 
tic. Results show that the actuator can be modeled as a proportional transfer 
element at up to 300 Hz (cutoff frequency). These results could be achieved 
only by using special guided roller bearings to reduce friction in the guides 
of the connecting rod. Magnetic bearings and electromagnetic actuators 
require a relatively large amount of space In relation to the magnitude or 
attainable forces; this drawback may be avoided by the use of hydraulically 
controlled chambers. 

Active chamber system. - The compact active chamber system, shown l^n 
figure 5, Is capable of generating very large forces and can thereby Influ- 
ence even large turbines weighing several tons (see ref. 15). The actuator 
device consists mainly of four cylindrical chambers arranged In a circle 
around the outer bearing housing. Each chamber Is sealed at the top and bot- 
tom by an elastic membrane. In order to decouple the forces In the x- and 
y-dlrectlons, the bearing housing Is supported against the membrane system by 
linear roller guides. The Influence of friction Is thereby reduced as we 1 . 
Alternative methods of support result from design variations. One variation 
supports the bearing via elastic rods - another, simply by chambers designed 
In such a way that they have an adequate transverse compliance. 
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Actuator pistons, as an alternative to membranes, have the disadvantage 
of possessing relatively large moving masses for the desired small displace- 
ments (less than 0.02 In., or 1/2 mm); this has an unfavorable effect on the 
transfer behavior at high frequencies. In addition, sealing problems occur, 
and friction forces appear on contact surfaces. 

From reference 16 a relation can be derived for the force generated as a 
function of the fluid pressure P and the radius of the membranes R: 



The factor Fg Isa function of the ratio of outer radius R to Inner radius 
r and the reference magnitudes P 0 and R 0 . Setting R/r - 3, P 0 » 1 bar, 
and R© - 1 cm, we obtain Fg ■ 35 N. If we then choose R - 10 cm and P = 10 
bar, the force generated would be F * 35 kN. The membrane system has to be 
designed In such a manner that both the stress caused by the pressure, and 
the displacements are In a reasonable range. Here, a wide range of variation 
is given by a specific choice of parameters. 

For Industrial applications the frequency range (more than 200 Hz) Is very 
Important. To Include the actuator system Into the overall system, with the 
target being to design an efficient control, the Input-output characteristic 
of the system Is necessary; It Is given In the x-dlrectlon as 


T(jco) 


F x (t) 

u x (t) 


(3) 


where u x (t) is the control voltage (Input) to the servo valve and F x (t) is 
the resultant force acting on the bearing housing. In the Laplace domain we 
can obtain the resultant force as a function of the system parameters by 
applying the continuity and the Bernoulli equations (under consideration of 
hydraulic losses), the dynamic behavior of the servo valve, and the force equi- 
librium at the membrane, assuming neglect of the oil compressibility. The 
resulting force equation Is given by 


where 



or by transformation Into the time domain the force equation becomes 


(4) 

(5) 

( 6 ) 
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(7) 


F * f V “ f 2 " f V A C T X B 


A C V X o 


A*C r 



C M X B 


where fy can be obtained by the differential equation of the servo valve 
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a),, 



A*K 


V 


pq 


( 8 ) 


or 


m v f V 


* d v f V 


* c v f v ■ 


b U 

V V 


(9) 


The valve constants K y, K pq , the valve elgenfrequency u v , and the valve damp- 
ing coefficient E v are constants normally given by the manufacturer of the 
valve; A* is the active membrane area, Is the membrane stiffness, Cj is 
a constant describing the oil Inertia, Cy and Cp are constants describing 
fluid losses, x 0 Is a characteristic flow velocity, U v Is the control voltage 
to the valve, and x B Is the deflection of the bearing In the force direction. 


Equation (7) describes the force that Is applied to the bearing, and equa- 
tion (9) represents the dynamic behavior of the valve that Is stimulated by 
the control Input U v . 


An optimal design of a control system demands an adequate model of the 
complete system. The effectiveness of the actuator system will depend on the 
observability (which, in turn, depends on the measuring information) and the 
controllability of the natural vibrations that are to be Influenced (which 
depends on the actuators used and their locations). In order to check on 
these system properties and to design an effective control system, adequate 
modeling is necessary. 


MOOELING OF ACTIVELY CONTROLLED ROTOR SYSTEMS 
Modeling of Rotor Systems to be Controlled 

Modeling the rotor system as a hybrid multibody system (HMBS) Is a very 
efficient method of describing It. This has proven especially useful when 
active components are included. The HMBS model contains rigid bodies (e.g., 
bearing units) and elastic subsystems (e.g., rotors and blades), as determined 
by special needs (see refs. 11 and 17). The coupling of subsystems is accom- 
plished by special elements that are characterized by the respective force 
laws. Modeling the system by this procedure permits the consideration of the 
disk mass, the shaft mass, and the gyroscopic effects in a simple manner. It 
also allows a modular construction of the system equations without any loss of 
physical Insight. The procedure leads to a computer-oriented representation 
In which control concepts can be adequately considered. It has the advantages 
of low system order, versatility, and simple system adaptation (e.g., altering 
location of control forces acting on a flexible shaft). The equations of 
motion of such HMBS are best set up on a computer. A very efficient method Is 
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the direct evaluation of the principle of d'Alembert, which can be expressed 
for a system consisting of k substructures as 


r dm - df e ^ Sr ^ - 0 (10) 

' S 

where the Index 1 indicates the substructure, r Is the position vector to 
the mass element dm, r Is the acceleration vector of the mass element, and 
df e Is the force vector acting on the mass element. The symbol S Indicates 
the variation of the position vector rj (no variation of the time t!). 

Rotor systems normally consist of elastic shafts and rigid bodies such as bear- 
ing housings or foundations. Taking Into account the assumptions of a 
Bernoulli beam (cross sections of an elastic shaft remain planar), we can Inte- 
grate over the mass elements In the radial direction and end up with a beam 
element (disk element) as shown In figure 6. (In figure 6 three different 
types of coordinate systems that are used In describing rotating machinery are 
also Indicated.) As a result of this integration, equation (10) becomes 




Here the first term describes the virtual work of translation (In parenthesis 
Is Newton's second law of motion), and the second term considers the virtual 
work of rotation (In parenthesis, Euler's equation); Sr Is the vector of the 
small virtual displacement, and S<p Is the vector of the virtual rotation; 
a s is the vector of acceleration of the center of mass; w Is the vector of 
angular velocity; « Is the vector of angular acceleration of the beam element; 
dl Is the tensor of moment of Inertia; and 



is the skew-symmetric tensor of the angular velocities given In the system of 
coordinates (see fig. 6). Assuming rotational symmetry relative to the shaft 
axis, dl for the disk element can be expressed as 


dl = 


0 

0 



(13) 


where I - I »^*r?d 2 - §?-d 4 (z>dz. The torque vector acting on the beam 
" x y 2 16 64 

element Is represented by dMI . In the case of a rigid body, the Integration 
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can be carried out over the entire body, and equation (11) would be expressed 
In the same manner but the differential quantities would become m, I, f e , 
and M e . 


If we now Introduce an f x 1 vector q of the generalized minimal coor- 
dinates, where f Is the number of degrees of freedom, the position vector 
r and the vector of the angle of twist can each be expressed as a function 
of q as follows: 

r = r(q) (14) 

which leads to 

= IS sq 

and 


<P » ?(q) 


(15) 


which leads to 


S 9 - f^Sq 

Their time derivatives can be written 


which leads to 


and 


r 



9r 

at 


8r 3r 
9q “ 3< l 


(ft) 


8q q + at 


which leads to 


(16) 


(17) 


9m 8m 
3q " 3q 

From the relations in equations (14) to (17), equation (11) can be expressed as 



where 8r78q - 8r/3q Is the jcoblan matrix of translation and 8<p/8q - 
9«/3q Is the Jacobian matrix )f rotation. The Jacobian matrices serve as dis- 
tribution matrices and can be 'nterpreted as transformation matrices of the 
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forces and the torques acting on the system In such a way that the forces and 
torques are projected onto the directions of the generalized coordinates. He 
will see the usefulness of these matrices later. 

A description of the system as an HMBS demands the subdivision of the gen- 
eralized minimal coordinates Into two types of hybrid coordinates: those with 

only time-dependent degrees of freedom (e.g., motion of the rigid bearing 
units) and those with degrees of freedom that represent elastic deformations 
(e.g., bending vibrations of the shafts). In the case of a shaft, the latter 
type describes, by distributed coordinates, the motion of the beam element 
shown In figure 6. These coordinates are a function of time t and the actual 
location of the beam element covered by the coordinate z along the rotor 
axis. To separate these two dependencies, we apply the Rltz product method 

T u T (z) 0 fi 

r.(z,t) - Cu(z,t), v(z,t)] \ - T q,(t) q ,6 R (19) 

1 1 0 v ( z) 1 1 1 


The vectors u(z) and v(z) are admissible shape functions, with fj as the 
number of shape functions considered In describing the deformation of the 
shaft (1^ substructure). These vectors can be looked on as expansions of 
real deformations that are described by these functions. The closer these 
functions are to the eigenfunctions of the real system, the better Is the con- 
vergence of the solution. (Note that the eigenfunctions change with the rotor 
frequency.) Good results have been achieved by using eigenfunctions of the 
nonrotating rotor, which allow a low system order even at high rotor frequen- 
cies (depending on the gyroscopic Influence). The determination of these func- 
tions Is carried out In a separate calculation, for example, by using cubic 
spline functions, Hermlte polynomials, or applying finite element or experimen- 
tal methods. The number of elastic degrees of freedom fp that are chosen 
depends on the elgenfrequencies of the system Itself and the expected excita- 
tion frequencies. The equation of motion of the separate modules (substruc- 
tures) can be given In the usual form of mechanical systems 



The coordinates that are necessary for the description of the separate modules 
are Included In the vector of generalized coordinates . The force vector 


h i ■ Ell j Iv <z - v dz = L jT<z 3 >f 3 <21 

i 0 i 

Is the sum of all forces acting on the 1 th substructure. The forces are not 
only Imposed from outside the system, but they can also be caused by coupling 
effects between the subsystems. The submatrices appearing In equation (20) 
are displayed In table I for an axlsymmetrlc shaft Cu(z) - v(z>], which holds 
true with very few exceptions. 
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TABLE I. - SUBMATRICES OF ELASTIC SUBSTRUCTURES WITH RESPECT TO INERTIAL SYSTEM 


Submatri cies 

In f 1 uence 

(M ). = (H ). =< 
xx i yy i j 

L f 

p[A(z)u(z)u (z) + Mz)u’(z)u’ (z)]dzs 

0 )i 

Inert i a 



'Vi - 'V. • "A 

( 

71 

EI (z)u"(z)u" T (z)dz\ 

o X > 

Material damping 

/ 

(G ), = (G ). =<2 
... | 

Q 

L ) 

pi (z)u’ (z)u' T (z)dz\ 

X / 

0 )i 

Gyroscopic effects 

(K ) . = (K ). = < 
xx i yy i \ 

L ) 

El (z)u"(z)u" T (z)dz \ 
x / 

0 h 

Elastic forces 

(N ). = (N ). = < 
xy l yx l \ 

r L t ) 

[2QpI (z)u' (z)u' T (z) + Qd.EI (z)u'(z)m" (z)]dzl 

X } x / 

0 21 

Spin and 

internal damping 


With equations (16) and (17), the Jacobian matrix In equation (21) can be 
expressed 


J 


J 


3r ;) (q 1 > 
3q 


( 22 ) 


where 


r^ = [u(q,Zj> . v(q,Zj)] T 


From equation (19), we obtain 




u T( z 3 ) 

0 


0 

* T( V 


(23) 


The vector of generalized coordinates of the overall system Is 

q . [qj.qj, . • -.qj] ( 

where q e R f , f Is the number of degrees of freedom of the overall system, 
and k 1 s the number of substructures (flexible and rigid parts). The 
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Jacobian matrices for the overall system relative to subsystem 
ten as 




can be wrlt- 


(25) 


The structure of these Jacobian matrices Is very simple; J, contains only 
zero matrices, with one exception: , which relates to th£ actual subsystem 

(ref. 11). 

By formulating the system In such a way, the substructures can be combined 
systematically to synthesize the overall system. The Jacobian matrices serve 
to connect the elements. In other words, they distribute the actions of the 
Interconnected subsystems In the matrices and thus describe the overall sys- 
tem. The set of equations for the whole system exhibits the structure of ordi- 
nary mechanical system equations 

Mq + (D + G)q + (K + N)q = h q G R f (26) 


By Introducing the state space vector 


x(t) - 



(26a) 


where x6R n and n ■ 2f, equation (26) becomes 

x - Ax + Bu 


(26b) 


where A G R n,n and B G R n,r ; f Is the number of degrees of freedom; r, the 
number of actuators; and u, the control vector. 

Using this method to model systems consisting of rigid and elastic subsys- 
tems has proven to be very suitable, especially when considering active ele- 
ments. The advantage of this method, compared to others. Is the low system 
order resulting from always using a Ritz approximation, which takes the mode 
shapes of the nonrotating rotor as admissible shape functions. Using these 
functions for the Ritz series provided very fast convergence, even at high 
rotor frequencies. For this reason the calculation expense is always low, 
even for very complicated rotor structures. 


Modeling of the Complete System 

To design a control system. It Is necessary to Include the actuators In 
the system model. The differential equation must be represented In such a 
form that the right side of the equation contains only terms that are Influ- 
enced by the controller (disturbance forces can always be added easily). 

The differential equation, especially the resulting changes In the system 
matrices, must be tailored to the electromagnetic actuators acting both 
directly on the rotor and via the bearing housings, and also to the active 
chamber system. 
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Inclusion of electromagnetic actuators . - We start with the following 
differential equation: 

Mq(t) + Pq(t) + Qq(t) = Z\ (q,q,q,Q,T, t) (26) 

where Q Is the rotor frequency, T Is the oil temperature when using journal 
bearings, and t Indicates the time. By using electromagnetic actuators to 
apply forces directly on the rotor , 


h 


J 



6(z - z^)dz = J^(Zj)fj 


(27) 


where Jr Is the Jacobian matrix of translation for the subsystem flexible 
rotor ; fj Is the discrete force vector (here control force) acting on the 
rotor at the location zj ; and S(z - Zj) Is the Dirac function that permits 
the consideration of the discrete location of Zj (screening characteristic of 
the Dirac function). 

Alternatively, we can apply forces that act Indirectly (via the rigid 
bearing housings) on the rotor 


h B = J B f B 


(28) 


where Jg is the Jacobian matrix of the subsystem rigid bearing element and 
fg is the force vector acting on this bearing element. Combining equations 
(1) and (28), we can formulate 


J 


T 

B 




O' 

k. 


x 

.y. 


+ J 



j t - 

Vb 1 


kJnJ„q + k 3 j;i 


A 


(29) 


In this formulation we took advantage of equations (14) and (16) to express 
the x-, y-movement of the bearing by using the vector q of the generalized 
coordinates and the Jacobian matrix of translation belonging to the subsystem 
bearing element, which is expressed as 


r 


B 


C*,y]J - J B q 


(30) 


When a magnetic bearing acting directly on the rotor is applied at location 
zj, equation (30) becomes 

r, = [u(Zj ,t) ,v(Zj ,t)] T - J R<Zj)q <31) 

The equation of the complete system can be achieved by Integrating the force 
equation (29) Into the differential equation (26a) 

Mq + P(Q)q + [Q(fi) - kjJj]q - Bu (32) 
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where B - kjknJR * s the control matrix with kq as the gain coefficient of 
the power amplifier, and u Is the control vector containing the Input 

voltages to the power amplifier. The negative sign of the term kjljg 
Indicates the negative stiffness produced by the magnetic actuator. 

Inclusion of active chamber system . - The active chamber system (hydraulic 
actuator) Is described by the force equation (8). We Introduce the 
abbreviation 


f - f v - mx' B - dx B - CgXg (33) 

and use the equation of the servo valve described by 


(9) 


as already given. The system order Increases by one degree of freedom for 
each servo valve, so the vector of generalized coordinates, q, has to be 
enlarged by the two additional coordinates, fy x and fyy (for the servo 
valves serving the x- and y-dlrectlons respectively): 

q - Cq T . fy x . fvy3 (34) 

The resulting differential equation, In conjunction with equations (31), (33), 
and (9), can be expressed as 


o * 1 

a o 
< 

o o , 

q + 

POO 
0 d v 0 

q + 

Q 0 0 

0 c v 0 

1° 0 m vJ 


O 

o 

Cl 

< 

i 


o 

o 

o 

< 


' f v, 

-Sx B 

-dx B 

- C B X B 

+ bjjj 

' U vx' 

1 — 

-h 

< 

-Sy B 


“ C B y B. 

V V 

U vy. 


(35) 


The coefficients m, d, and Cg describe the additional mass, damping, and 

stiffness effects caused by the actuator system; Jg Is the enlarged Jacobian 

matrix of translation for the subsystem bearing, and J y Is the Jacobian 

matrix belonging to the subsystem servo valve. Those terms of equation (36) 
that are multiplied by J g are linearly dependent on the system coordinates 

and thus contribute to the homogeneous terms of the differential equations 
(dynamic behavior of the complete system). 

The radial movement of the bearing, described by xg and yg, and the 
servo valve dynamics, described by coordinates f y x and fyy, can be expressed 
by the enlarged Jacobian matrices of the subsystem (ref. 11) 
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( 36 ) 




and 


1 Vx 


L VyJ 




(37) 


The final formulation of the differential equation of the complete system can 
be given as 


[m + OgJgmjq + [ P + J B J B d ] q + b + J B J B C B " J B J v] q " Bu 


(38a) 


or 


Hq + Pq + Qq - Bu (38b) 

One can see that the coupling of the mechanical system, the fluid system, and 
the electrical system can be managed very elegantly. The Jacobian matrices 
that are Involved serve as the "distribution-matrices.” They distribute the 
Interaction of the Involved subsystem. 

The Input by the actuator Is characterized by the term Bu where B Is 
the control matrix, which depends to a high degree on the actuator locations 
chosen, and u Is the r x 1 control vector with r as the number of actua- 
tors Involved. 

If we use a linear combination of all Imaginable measurements (displace- 
ments q, velocities q, accelerations q. and Integrated displacements qp, 
the control vector can be expressed as 

u = -Kc[q[,q T ,q T ,q T ] =» -KCx 
or 

u - - [Kj ! Kp ! Kq ! Kfl] Cx 

where Kt Is the gain matrix relative to qj - Jqdt; ICp Is the gain matrix 
with respect to the displacement vector; Kq Is the gain with relation to the 
velocities; Ka Is the feedback matrix taking Into account the accelerations; 
and C Is the measurement matrix, which takes Into account that often not all 
displacements with their time derivatives and Integrals are known. The rela- 
tion between an enlarged vector 


(39a) 

(39b) 


x 



W.q 1 ] 


T 


and the measurement vector y can be given as 
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( 40 ) 


(41) 


y = Cx 

With equation <39a>, equation (38b) can be formulated as 

<M + BK A C)q + (P + BKpOq + (Q + BKpC)q + BKjCqj - 0 

For several reasons It Is more convenient to transform equation (42) Into the 
state space representation. By Introducing a state space vector 

x . [qJ,q T .q^] T <«> 

(where x Is different from the normally used state space vector because it 
has been enlarged by the vector qj, which contains the Integrated displacement 
coordinates) equation (42) becomes 


or 



' 0 

E 

0 

0 

0 

E 

-M^BKjC 

-M^ 1 (P + BKpC) 

-H a ] (Q + BKpC) 


(44) 



(45) 


Equation (42) shows that the mass matrix M A U a function of the acceleration 
feedback characterized by K A : 


M a - H + BK a C <46) 

In the case of acceleration feedback - this depends on the software available - 
using the second order equation (42) for optimization of the controller can be 
advantageous because of the controller-dependent mass matrix 

m a - m a oca> < 47 > 


Discussion of Several Cases 

ProDortlonal derivati ve (PD) controller . - In this case, where ICp w 0, 
Kn * 0, Ka - 0. and Kr - 0, M A Is constant; however, because Kt - 0 (no 
Integral term), the system can be reduced to the usual form (the fi x 1 vec- 
tor qr In the state space vector x can be canceled). This means the 
upper f i rows and the fj columns at the left side of the system matrix A 

can be canceled. 

Proportional derivative-acceleration (PDA) controller . - Here, Kp * 0, 
Kn * 6, Ka # 0, and Kt - 0. "The reduced order system (PD controller) Is 
still valid here, but the mass matrix Is a function of the acceleration feed- 
back [M A - M(Ka)]. 
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Prooortlonal Integral derivative-accel eration (PIDA) controller . In^ 
this case, all submatrices'" appearing In K (see eq. (39)) are n ° nze ™; 
for calculation (optimization and simulation), the enlarged system described 
by the differential equation (45) - full order system - has to be taken into 

account . 

From the complete description of the system, the design of suitable con- 
trol concepts can be achieved. The effectiveness of the controller, In connec- 
tion with the actuators and sensors used, depends to a high degree on the 
controllability and observability conditions of the system. These system prop- 
erties will be discussed In the following section. 


CONTROLLABILITY, OBSERVABILITY AND SPILLOVER 
Controllability and Observability 

Controllability and observability are two of the most Important factors 
In the theory of dynamic systems. They play a very Important role when des g 
ing a control system and give an Insight Into the physical problem. S mply 
out. controllability means the ability to adequately control the dynamic behav- 
ior , orThe _ s7Ttenrthrough the actuators used. Observability means that appro- 
priate sensors exist at appropriate locations such that the pertnent dynamic 
behavior can be detected. The controllability and the observability J h J h r 
system are mainly determined by the chosen actuators and sensors and their 
locations. Controllability and observability can change as a function of 
rotor speed, stiffness and damping of the bearings, oil te JPerature of the 
fluid film (In the case of journal bearings), and so on. Much software Is 
available to check on these system properties. For tlije-lnvari ant mechanical 
systems, the Hautus and the Kalman criteria are generally used (see refs. 17 

and 18). 

The controllability and/or observability can be studied via the mode 
shapes (eigenvectors) of the system. This can be explained on a s mple system 
whose mechanical model Is shown In figure 7. When one uses control forces, 
the controllability Is Indicated by the amplitudes of the various mode shapes 
at the actuator location. This means the larger the amplitude, the higher the 
controllability. For Illustration, the mode shapes up to the fourth order are 
displayed In figure 8 for the system shown In figure 7 for two different rotor 
frequencies (Q - 0 and Q - 100 Hz). Only one set of modes appears when 
q , o (fig. 8(b): symmetrical system, decoupling In x- and y direction). 

The modes split at rotor frequencies Q * 0 Into forward and backward modes 
(figs. 8(c) and (a) respectively). 

By applying forces that act directly on the rotor (with a magnetic bear- 
Ina) all natural mode shapes of the rotor could be Influenced sufficiently If 
the forces acted at locations where the amplitudes of the modes have their max- 
imum value Unfortunately, this maximum occurs at different places for dlffer- 
eTm^des and the modes c^n also be Influenced by the rotor frequency (compare 
", “d KcESSird modes at Q - 100 Hz, gyro influence). The optimal 

actuator location depends on one's objectives, and It will always be a 
comproml se. 

Now consider the rotor to be supported by roller or journal be J r1n 9*; . 

The vibration behavior Is to be Improved by control forces via bearing housings 
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(Indirectly acting on the rotor). In this case the effectiveness, that Is, 
the controllability, Is characterized to a large extent by the stiffness of 
the bearings and also by the bearing support stiffness. 

For some Idea of these Influences see figure 9; it shows a measure of the 
controllability for the first and the third forward modes of the system In fig- 
ure 7. As a measure of controllability, one can use the square or any other 
suitable norm of the bearing displacement at the actuator location. In each 
case the controllability Is plotted as a function of the bearing stiffness 
C|< and the bearing support stiffness c a relative to the stiffness c w of 
the rotor Itself. In figure 9(a) the controllability can be seen to decrease 
with reduced C|<. Physically, this result Is very understandable; the lower 
the stiffness between bearing and shaft, the smaller the effectiveness of con- 
trol forces acting via the bearings. The behavior with respect to the outer 
bearing support stiffness Is the opposite. Although not shown, these results 
are similar for the second forward mode but are not valid In general. The con- 
trollability of the third forward natural mode (bearing displacement) Is plot- 
ted In figure 9(b). This behavior Is different from the preceding one because 
of Increasing influence of inertia with rising frequency of oscillation. 

For the example under consideration, good controllability of the important 
natural mode shapes (up to the third order) In the current frequency range is 
guaranteed If the stiffnesses satisfy the conditions 

c k c a 

— > 1 and < 1 (48) 

c * c w 

Investigations of the observability are similar, and for the system under con- 
sideration the results are Identical. 


Spl I lover 

If elastic components are Involved In the system, an exact system descrip- 
tion requires an Infinite number of shape functions. In practice, only a 
finite number of modes can be considered, and In fact, only a few are needed 
for adequate control. The consequences of this Incomplete system description 
are spillover effects. 

Spillover effects can appear on two different levels. On the one level, 
the measurements con tai n both the modes considered In the model and also the 
modes that have been disregarded. Both sets of modes Influence the signals at 
the measurement locations; this leads to observation spillover . On the other 
level, control forces may destabilize modes of higher order that have not been 
accounted for In the control design, thereby leading to control spl 1 lover . 
Usually, Internal damping Is assumed to be large enough to prevent Instabili- 
ties due to spillover effects. Because this Is not true In every case, It Is 
prudent to consider these effects In the control design. Theoretical Investi- 
gations lead to the collocation condition (refs. 5 and 20) 

B - C T (49) 

where B Is the control matrix (eq. (26 b)) and C Is the measurement matrix 
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(eq. (41)). The condition given by equation (49) requires that actuators and 
sensors be located at the same place. It also requires only velocity feedback; 
the number of signals Involved In the feedback must be equal to the number of 
actuators used. 

For economic reasons, using more signals than actuators In order to 
receive as much Information as possible about the system seems quite reasona- 
ble. In fact, an Improvement of the dynamic behavior is often possible if we 
use more sensors than actuators. To avoid Instabilities, the neglected mode 
shapes In the system description must be neither observable nor controllable. 
These requirements lead to a practical criterion for the choice of the actua- 
tor and sensor locations, which Is expressed by the following equations: 

S | v ^ (z<-) | ♦ min 1 > n n <50) 


and 

£ | v . (z A ) | ■» min - i > n n <51 ) 

where z$ is the position of the sensors, z/\, the actuator location, and n n , 
the number of shape functions considered in the system description. Reference 
5 provides hints on managing the simulation of such a system. 

Minimal spillover can be expected if at the actuator locations the ampli- 
tudes of the higher modes tend to zero. Often, suitable positions may be 
found directly by evaluating the eigenfunctions (fig. 10). For systems with 
large mass concentrations, the high-order natural vibrations have almost no 
corresponding rigid body motions. Only the low-order relevant vibrations con- 
tain significant rigid body movements In addition to elastic deformations. 
Therefore, we recommend placing the actuators and sensors at the points of 
mass concentration where, In practice, the higher order eigenfunctions have 
none or only small displacements. When we take a closer look at the mode shape 
with a frequency of 2100 Hz (fig. 10), we can see that to avoid spillover, the 
sensor should pick up the bearing movement and not the movement of the rotor 
at the bearing location. 

Spillover effects may also be due to sensor locations and actuator loca- 
tions not being the same. For example, two rotor modes are Illustrated In fig- 
ure 10. For the first mode, the rotor moves In the same direction at both the 
sensor location and the actuator location. For the second mode, the rotor 
moves In opposite directions at the two locations. This means that, for first 
mode motion, a positive x — d 1 splacement at the sensor requires a negative 
x-force at the actuator to counteract it. However, for the second mode, a pos- 
itive x— d 1 splacement requires a positive x-force to counteract it. That is, 
the relation between sensor motion and required actuator force is opposite for 
the two modes. If the controller gain Is assumed constant at all frequencies, 
and the controller Is programmed to suppress first-mode vibration, then the 
second mode will be unstable. This occurs because the actuator response Is to 
Increase the second mode motion rather than to suppress It. 

Spillover can also be shown analytically as follows. Consider the case 
when velocities x$ and y$ are taken at location z$ with the gains 
k x ® ky = k and without coupling between the x- and y-di rectlons . Then for 
control forces generated by a magnetic bearing located at plane M, we obtain 
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(52) 



q 


X 


-kJ(z s > 



For a demonstration of spillover effects, only two modes will be considered 
(that Is, two In each direction); see equation (19). The number of degrees of 
freedom f Is 4. The vectors of the admissible shape functions become (com- 
pare wl th f Ig . 11) 


u(z) = v(z) 




. TfZ I 

V 


sin — 


* 

. 2irz 

, U 2. 


sin -j— 


(53) 


The sine functions In the sample considered are the eigenfunctions of the non- 
rotating rotor. For Inclusion of the force vector (eq. (21)) In the dif- 
ferential equations given by equation (20), the force vector can be expressed 


h 


-kJ T (z M )J(z $ ) 



(54) 


With equation (53), the Jacobian matrices appearing In equation (54) can be 
given explicitly as 




U M1 U M2 0 0 


0 0 U M1 U M2 


and 


u s , U S2 0 0 


0 0 U SI U S2, 


(55) 


If we take z^ = 1/4 and z$ - ( 3/4) L , we obtain u^i - u$i - yfT/2 = 
a > 0 and um 2 * -u$2 » 1 - b > 0. The spillover effects can be explained ana- 
lytically If we set the rotor frequency Q = 0. In this case, equation (20) 

Is decoupled in the x- and y-dlrectlons (see table I). Assuming the runner 
consists of a thin disk with mass m and radius r, the thickness b of the 
disk Is much less than the length L of the shaft, and the shaft is massless, 
then equation (20) can be written In the x-dlrectlon as 



1 0 


*1 


V 

-ab 


q i 


C 1 

0 


q l 

m 




+ k 

ab 

-b 2 


i« 

+ 

0 

r 


q 2 


o 

0 L 2 J 


q 2 



q 2 


V 

c 2 



(56) 


Note that equation (56) represents only half of the system. The other half 
(In the y-dlrectlon) Is the same because of rotational symmetry. 

The stability of the system can be verified by the characteristic equation 
relative to equation (56) 

AmX 4 + (ka 2 A - mkb 2 )\ 3 + (c 2 m + c 1 A)X 2 + (c £ ka 2 - c ] kb 2 )X + c,c 2 - 0 (57) 
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where A - 2ir 2 r 2 m/L 2 . By 
ble if the coefficient of 
of x is 


the Stodola criterion (ref. 18), the system Is unsta- 
X and/or X 3 Is less than zero. The coefflcent 


a] - k(a 2 C2 - b 2 q) 

and the coefflcent of X 3 Is 

33 = k(a 2 A — mb 2 ) 

From table I and equation 53, we obtain 

nL 


dz-©EI© 


and 


d 0 




c 2 - li 


6© EI s1n 2 ^dz.l6(0El(t) 


d 0 


(58a) 


(58b) 


(59a) 


(59b) 


(60) 


(61) 


Combining equations (59a) and (59b), we obtain 

C 2 - 1 6c i > 0 

Since b 2 * 2a 2 >0, It follows that If k > 0 

ai > 0 

As a result of equation (61), we know, for positive k, that the system can 
only be unstable If a 3 < 0. If L/r > *, this condition Is fulfl led and the 
system becomes unstable. (If k < 0, then aj < 0 and the system Is unstable 
for any value of L/r.) 

CONTROL CONCEPTS AND CONTROLLER OPTIMIZATION 

The design of linear and nonlinear controllers, both analog and digital, 
has been treated in great detail In the literature. There are numerous methods 
and concepts for the design and realization of controllers. The choice of a 
concept is determined mainly by the aim and the physical real 1 ties of the plant 
(l.e., the open loop). In the following, only a time Invariant and linear sys- 
tem is assumed: 

x(t) - A(Q,fi,t,T)x(t) + Bu<t) A € R n,n and B 6 R n * r 


and 


y(t) - Cx(t) 


C G R m » n 


(62) 


Here A Is the system matrix; B, the control matrix (dependent on the actua 
tor locations and actuator type); C, the measurement matrix (dependent on t e 
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measurement locations and the type of sensors used); n - 2f, the system order, 
where f Is the number of degrees of freedom; r, the number of Independent 
forces generated by the actuators; m, the number of signals used to realize 
the feedback (measurement coordinates); Q, the rotor frequency; ft, Its time 
derivative; and T, the oil temperature. 

Depending on the Information used, two concepts may be distinguished: 
state feedback 


U< t ) - K S (Q,T)x(t) K S € R r » n (63) 

and output feedback 

u(t) = KQ(Q,T)y M (t) Kq G R r . m (64) 

This means that the control vector u Isa linear function of either the sys- 
tem state vector x(t> or of the system output, that Is, the measurement vec- 
tor yM<t), which Is directly available from the system sensors. 


State Feedback Control 

In designing a state feedback controller, the following methods are gener- 
al ly used: 

(1) Optimization according to the quadratic Integral criterion (or per- 
formance Index): 


J * f (x^Qx + u T Ru)dt ■* min (65) 

J 0 

This criterion takes Into account all coordinates for describing the sys- 
tem represented by x and for the system output given by u as well. By 
using this approach, the design problem Is reduced to the task of solving the 
algebraic Rlccatl equation (ref. 11). The weak point In this procedure Is the 
adaption of the weighting matrices Q and R for the specific problem. 

(2) Choice of the eigenvalues Xj of the closed control loop (pole 
assignment) : 

CXi E - (A - BK)]q - 0 (66) 

In contrast to the previous method, the choice of suitable eigenvalues 
(poles) Is the problem with this approach. With the selection of eigenvalues, 
a specific system behavior Is enforced, so In the case of poorly adapted val- 
ues, extremely high control forces may appear. The Influence of gyroscopic 
effects may even worsen the system dynamics. Hints on choosing suitable poles 
are given In reference 21. 

(3) Modal state control 

This method allows for a shift of one or more eigenvalues. This can be 
especially useful for rotor systems that, for whatever reasons, are run In the 
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vicinity of resonances. With the control concept given In the following, how- 
ever, one can only shift as many poles as there are actuators (refs. 5 and 21). 
In the control design state, the system (eq. (62)) Is transformed Into 

x(t) = Ax(t) + HBuU) = [A - K|m|]x< t) (67) 

The matrix A Is real and has a diagonal block structure (ref. 11). The gain 
matrix appearing In equation (63), which Is needed to realize results, Is 

given by 

K S - B+TKmT" 1 <68) 

with 


B + - (bTbHb 1 <69) 

being the pseudoinverse of the column regular control matrix B from 
equation (62). 

(4) Combined state feedback control 

This kind of control Is useful in cases where the application of a con- 
troller that was designed according to method (1) or (2) does not yield suffi- 
cient damping of some natural vibrations. Here the additional application of 
a modal feedback control according to method (3) can Improve the system behav- 
ior significantly. The additional effort for this case Is negligible since a 
summation of gain matrices affects only the gain coefficients, whereas the 
structure of the controller remains the same. 


Output Feedback Control 

With a limited number of measurements, the possible Influence on system 
behavior Is restricted. However, the obtainable results are still sufficient 
for most applications. Significant Improvement by means of state feedback 
often requires a large number of sensors, which In most cases Is not practical. 

The design of an output feedback control always Implies a parameter opti- 
mization, that Is, an optimal tuning of the coefficients of the feedback matrix 
In equation (64). If we design an output feedback with the quadratic quality 
criterion, the quality functional can be acquired by solving the Ljapunov 
matrix equation. Utilizing all symmetry characteristics reduces the number of 
variable parameters significantly, thereby decreasing the calculation expense 
as well. 

It Is Important to note that In contrast to the Rlccatl controller, dif- 
ferent Initial conditions x(t 0 ) - Xq will result In different optimal con- 
trol matrices Kq. The Initial conditions may be set explicitly In such a way 
as to enforce consideration of some critical natural mode shapes (e.g., by 
equating x 0 to a specific eigenvector characterizing a natural mode shape). 

Under certain circumstances (l.e., when fully observable and controlla- 
ble) the output feedback control allows for shifting of specific eigenvalues 
as well. If the modal description of the system Is as In equation (67), the 
feedback matrix Kq In equation (64) becomes 
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wl th 


K o - b+tkmHc- 


(70) 


C+ = C T (CC T H (71) 

being the pseudoinverse of the row regular measurement matrix C in equation 
(62). In general, the matrix C Is not square so it cannot be inverted; thus 
the pseudoinverse indicated In equation (71) must be used. This means that 
the other poles are shifted as well, and a check on the stability and system 
behavior after the shift of the poles is required (refs. 21 and 22). 


EXAMPLES AND RESULTS 

Knowing how to construct a controller and which type of controller to 
chose depends to a high degree on one's objectives. In this section, several 
examples will be given to demonstrate the design procedures. 


Magnetic Suspension of an Epitaxy Centrifuge 

The information in this section Is taken from reference 1. The rotor 
system was developed for application in liquid-phase epitaxial growth of very 
thin semiconductor layers. In order to obtain layers of high quality, a very 
smooth rotation is required in a reactor that is leakproof even under ultra- 
high-vacuum conditions. An active suspension of the rotor was necessary in 
order to absolutely exclude contamination by lubricants or by wear. 

The resulting epitaxy centrifuge, supported without contact, Is shown In 
figure 12. The upper part of the apparatus contains the rotating crulclble 
where the epitaxial layers grow. The crucible can be heated up to 1000 K by a 
furnace. The supporting unit Is concentrated In the lower framework. The 
rotor spins without contact In the vacuum tube. The bearings, the attitude 
sensors, and the stator of the electric drive have to be outside the tube to 
prevent contamination. Bearing forces, attitude signals, and driving torque 
are transmitted through the walls of the tube. Two emergency bearings prevent 
damage of the rotor In case of a magnetic bearing failure. 

Equations of motion and state . - The mechanical model that Is the basis 
of the system described Is shown In figure 13. The elongated rotor Is assumed 
to be rigid. Asymmetries of the rotor are restricted to small dynamic and 
static unbalances only. The driving torque acts about the rotor axis. The 
radial bearings exert discrete horizontal control forces; the axial bearing 
force acts along the rotor axis. The axial bearing's vertical force component 
compensates for the rotor weight m»g;, Its horizontal components act like an 
elastic restoring force. The deviation of the rotor from Its vertical refer- 
ence position Is described by the given vector r(t) and the small Inclina- 
tion angles a and 3, as Indicated In figure 13. Linearization of the 
equations of motion Is justifiable. As a consequence, the horizontal motion 
Is decoupled from the vertical motion. 

Because the rotor Is assumed to be rigid, the four degrees of freedom of 
the radial motions can be described by the four coordinates x and y 
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(translational) and a and ft (rotational). The vector of generalized coordi- 
nates for describing the dynamic behavior of the rotor can be taken as 

q - [x.ft.y.-*]T ( 72) 

where x and y are radial displacements of the rotor at point S shown in 
figure 13. The negative sign of the angle a can be explained by the dis- 
placements In the y-d 1 rectlon that are caused as a result of angular movement 
a. These displacements at the locations Z] = L u or Z 2 = L 0 are positive 
If a Is negative. 

From equation (19) we obtain the vector of the admissible shape function, 

u(z) » v(z) » [1 ,z]T (73) 

or the displacement vector as a function of the coordinate z along the rotor 
axis , 


r(z,t) - \- - jji? - -1 Cx,ft,y,-a] T 
LO CM zJ 


(74) 


or In connection with equation (14), 

r(z,t) 


8r 

3q 


J T q 


(75) 


with Jj as the Jacobian matrix of translation. The Jacobian matrix of rota- 
tion appearing In equation (18) can be expressed 



"3a_ 

3a 

3a 

3a 

8ci> _ 8[a,ft]^ 

8x 

3ft 

3y 

3(-i) 

8q 8q 

8ft_ 

M. 

ai 

jL- 


-3x 

3ft 

3y 

3(-a)- 


= J 


(76) 


Radial forces acting on the rotor . - These are forces generated by the 
active elements and the weight of the rotor at the locations z ■ Ly . S , Lq , L/\ 
(fig. 13). The force vector of equation (21) can be given as 


h 


J T T(L u ,f Mu 


+ J T <L o >f Mo * J T (L A )f A * J r‘c 


(77) 


and f u and f u are given by the magnetic bearing equation 

MU "0 


Mu ,0 


L f yJ 


u.o 


0 


syj 


U.O 


Mj 


1. 


L^J 


U.o 


(78) 


or 

f. *> k^Jj(Zj)q + k^lj 
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The radial component of the axial magnetic bearing can be considered as a 
restoring force 


■ -wv (79> 


with k m as a bearing constant. 

Eventually, the torque t c produced by the couple of forces consisting of 
weight G and the axial force Fm for compensating the weight can be given 
as 



t c - -mg 


a(L A - S) 


3(L a - s> 


-mg(L A - s)J R q 


(80) 


Combining equation (78) with equation (80), we obtain 

" - [W'VVV * J ? a o )J T a o) k s - VJVW - n 9 <L A - * »“ 


(81) 


where 


u = t i xo» ^ xu * ^yo* ^yu^ (82) 

is the control vector built up by the control currents In the given order. 

With regard to equation (78), the control matrix B has to be 
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The matrices of the system 

Mq + Pq + Kq = Bu (84) 

can be calculated from table I and the vector of the admissible 
shape functions (eq. (73)) as follows: 
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Introducing the state vector x = [q ,q ] , we obtain the state equation 


x - Ax + Bu 


-M -1 K -M -1 P. 


and 


M -1 B. 


( 85 ) 


( 86 ) 


The essential values of the design parameters of the centrifuge are given In 
table II. 


TABLE II. - DESIGN PARAMETERS OF CENTRIFUGE 


Coordinates, m 

Lo 

L U 

s 

l a , 

Moments of Inertia, kgnr 

Transverse, Iy 

Polar, I 2 

Mass , m, kg .... 

Magnetic bearing coefficients, 

ki, N/A 

k s , N/m 

k m , N/m 

Gain of power amplifiers, kQ, A/V 


0.314 

0.126 

0.458 

0.448 

0.060 

0.00421 

0.458 

13.0 
. 6150 
. 0.1 
. 0.25 


If we measure all coordinates represented by the vector q and differen- 
tiate all the signals, we know the whole state vector x. In this case, the 
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system Is fully controllable and observable. Using a Rlccatl controller to 
optimize the controller with regard to the Integral criterion given by equa- 
tion (66), we obtain the damping represented by the real parts of the eigen- 
values of the system shown In figure 14. 

The weighting matrices for the results displayed by the dashed lines were 
chosen as 


Q = 



O' 

0 . 



4 4 7 

E £ q Q - 2.5x10' 


(87) 


For these results the rotor frequency chosen to design the controller was 
0 = 150 Hz. Figure 14 shows that the rotor Is unstable up to a rotor fre- 
quency of 30 Hz, since there Is negative damping associated with one of the 
eigenvalues below this frequency. To avoid Instability, we can use an adapt- 
ive controller. This means that the controller always has to be adapted to the 
actual rotor frequency. The attainable results are displayed In figure 14 by 
the solid lines. 


Control of an Elastic Rotor 

A rotor consisting of an elastic shaft with a rigid body at the top sup- 
ported by ball bearings was Investigated (fig. 15). The shaft diameter varies 
with axial coordinate z; the cross section of the shaft Is constant by sec- 
tors. To Improve the dynamic behavior of the rotor system, a magnetic bearing 
was used (fig. 15). 

The essential values of most design parameters are given In table III. 


TABLE III. - DESIGN PARAMETERS OF ELASTIC 
ROTOR SYSTEM 
[See fig. 14.] 


Coordinates, m 



s . 

Moment of Inertia of rigid body at end, 


0.0425 

0.2565 

0.5 

0.0652 

kgm 2 


Transverse, I^y .... 

Polar, I r z 

Mass of rigid body, mp, kg 
Bearing stiffness, N/m 

c U 

Young's modulus, E, N/m 2 
Mass density, p, kg/m 3 . 


0.0986 
. 0.131 

. 14.42 


2 .OxlO 7 
1 .0x10® 
2.06xl0 1] 
7.85xl0 3 


The equation of motion can be achieved by using the theory given In the 
section Modeling of Actively Controlled Rotor Systems. 
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In the first step, the mode shape functions have to be determined. The 
determination of these functions Is carried out separately; for example, by 
use of cubic spline functions, Hermlte polynomials, finite elements, or experi- 
ments. The results of this determination are shown In figures 16 and 17. Fig- 
ure 16 displays the elgenfrequencles with respect to rotor speed, up to the 
fourth order. Note that all the curves split Into one forward and one back^ 
ward vibration mode. 

Because only the forward vibration modes can be excited by Imbalances (In 
most cases, the dominant excitation In rotating machinery), these modes must 
be considered relative to the rotor speed. For this reason, the runup line 
(» * Q) Is plotted In figure 16. Each Intersection between the runup line and 
a curve of the forward elgenfrequency signifies a critical speed of the rotor. 
In the plot shown, two critical speeds can be recognized: o>i - 12 Hz and 
gj? - 330 Hz. In order to see the effect of the chosen actuator or sensor loca- 
tion on the controllability and observability of the system, the mode shapes, 
up to the third order for three different rotor frequencies (fl - 0, 150, and 
300 Hz, are plotted In figure 17. The actuator location Is marked by M. The 
mode shapes and their frequencies are obviously functions of the rotor speed 
(gyroscopic Influence). A close look shows that the controllability of the 
third forward and the third backward modes at Q* 150 Hz and the third back- 
ward mode at Q = 300 Hz Is tending to zero; this Is evidenced by the small 
amplitude at location M. If these modes are to be Influenced effectively, 
the actuator location must be changed. The results attainable by state and 
output feedback control will be discussed In the following section. 

State feedback control (Riccatl controller) . - The control design is 
based on the quadratic Integral criterion shown In equation (65): 

J = j J (x T Qx + u T Ru)dt -► min (88) 


This leads to 


j . ■ i xl Px„ 
opt 2 o o 

where P is the solution of the algebraic Riccatl equation 

A T P + pa - PBR"1 B t P +0=0 


(89) 


which. In connection with the criterion of equation (88), always supplies an 
optimal controller (gain matrix K) 


u = -Kx 


R-^P 


(90) 


A significant amount of software Is already available to solve the algebraic 
Riccatl equation. 


With the vector of generalized coordinates 

i T 


k T - 


(91) 
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and the state vector 


x = Cy T ,y T ] 

the weighting matrices are chosen as 


Q 


0 0 0 

0 0 0 

o 0 Q u 

0 0 0 


0 ■ 
0 
0 


and R 


Q 


vj 


E 


Only the velocities are taken Into account. If we take 

Q u - Qv - d1ag{l .4xl0 7 , 9xl0 6 , 6xl0 6 } 
we obtain, for a rotor frequency of 150 Hz, the gain matrix 


(92) 


(93) 


(94) 


-0.43xl0 4 

-0.35xl0 5 

0.40xl0 6 

-0.81xl0 4 

0.29xl0 6 

0.27xl0 6 


k , 

k 3 

0.81xl0 4 

-0.29x10 6 

-0.27xl0 6 

-0. 43xl0 4 
-0. 35xl0 5 
0.40xl0 6 


~ k 3 

k l 

0.28xl0 4 

0.33xl0 4 

0.84xl0 3 

0.43xl0 3 

0.28xl0 3 

-0.93x10 2 


k 5 

k 7 

-0.42xl0 3 
-0.28xl0 3 
. 0.93xl0 2 

0.28xl0 4 
0.33xl0 4 
0.84xl0 3 . 


- k 7 

k 5 


(95) 


Because the rotor Is symmetric, the gain matrix can be represented 



(96) 


It is striking that the gain matrix applies not only to -(ks), but also to 
artificial gyroscopic (ky) and restraint (k]) forces as well as to nonconserva- 
tive forces <k 3 ) . The damping equivalents with respect to rotor frequency are 
plotted In figure 18. 

In this figure one damping curve (dashed line) Is negative (the corre- 
sponding vibration mode is unstable!). This Is caused by spillover effects . 
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To design the controller, only three mode shapes have to be taken Into account; 
however for the model to simulate the system equations, four mode shapes were 
considered. Note that the design criterion guarantees asymptotic stability 
only for the model that was the basis for designing the controller. Figure 18 
also shows that one of the damping curves approaches zero near Q = 350 Hz. 

This lack of damping Is due to the controllability condition of this mode 
shape at this rotor frequency (fig. 17). We may conclude that the state feed- 
back controller as designed Is useless. A reasonable state controller would 
be much more expensive, so an output feedback controller will be discussed as 
an alternative. 

Output feedback control . - In contrast to the optimization procedure of 
the Rlccatl controller leading to an optimal controller, the structure of the 
output controller will be fixed at the outset. For the feedback, only directly 
measured state coordinates are expressed by equation (62). By taking Into 
account this measuring equation, the control vector can be written 



u * — Ky = — KCx * — (K-jl^) 


10 



(97) 


As a criterion for design of the controller, the quadratic Integral criterion 
will be used: 


which leads to 



opt 


1 A* 

2 o o 


(98) 


where P Isa solution of A^P + PA - Q = 0 with A = (A - BKC) . 

The optimization can be carried out with a presupposed starting matrix 
Ko only If the system Is asymptotically stable. Now, by solving the Ljapunov 
equation (eq. (98)) and varying the gain coefficient appearing In K, the value 
of the criterion can be minimized- One weak point In this procedure Is that 
an Initial state vector x 0 is necessary. On the other hand, the procedure 
can be used for a higher valuation of a special mode shape (here, e.g., the 
first forward mode). This means that for x 0 , the eigenvector (mode shape) 
with respect to this vibration mode has to be determined. 

To reduce the calculation expense, we can take advantage of the special 
structure of K (properties of symmetry) 


K - (K^) 



(99) 


When compared to the results of the Rlccatl controller (eq. (95)), equation 
(99) shows that to achieve an optimal controller, forces that act like gyro- 
scopic forces (ky) and nonconservative forces (k 3 > have to be added. 
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The gain coefficients obtained when using the same weighting matrices as 
given by equation (93) are 


kj = -1 . 9x1 o 5 
k 3 = 2.3xl0 4 
k 5 = 2.4xl0 3 
k 7 = -9.0xl0 2 


(100) 


The results, plotted in figure 19, show that In contrast to the Riccati con- 
troller, the system is stable at all rotor frequencies. The slight damping of 
one of the vibration modes Is still there, but this problem can be solved by 
changing the actuator location (changing the controllability). 


CONCLUDING REMARKS 

Active vibration control of rotating machinery Is being given more and 
more attention by research Institutes as well as by industry. The purpose of 
this report Is to present the problems confronted when applying active control 
techniques to rotating machinery dynamics. The success or failure of active 
measures Is determined by the availability of appropriate actuators; by model- 
ing of the entire system, including all active elements Involved; by position- 
ing of actuators and sensors; and by the control concept used. All of these 
topics have been addressed, and their special problems have been discussed in 
detail In this report. A survey of existing actuators as well as those that 
are still in the design stage is Included. 

A very efficient method - called "hybrid multibody systems" - was 
described and used to analyze rotor systems consisting of rigid and elastic 
subsystems. This method allows a modular construction of the system, which is 
very easy to handle on a computer and Is both systematic and clear. Further- 
more, control aspects can be adequately considered (e.g., simple system 
adaptions, with respect to actuator and sensor locations, or optimization 
strategies for designing the controller). 

Such Important aspects as controllability, observability, and spillover 
were discussed. A method of checking on these system properties was outlined, 
and examples were displayed. Then, the most frequently used control concepts 
were introduced and their strengths and weaknesses pointed out. Real applica- 
tions served as examples to demonstrate how to design an optimal controller. 
These examples Indicate a possibility for Improving rotating machinery by 
applying active vibration control. 
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APPENDIX - SYMBOLS 


A 

A 

A* 

B 

b 

b v 

C 

c 

c M 

Cr 

Ct 

c v 

c 

C V 

D 

d 

d v 

E 

F 

FV(S) ,F2(S) 
f 

f e 

df e 

fv 

G 


system matrix 

moment of inertia 

characteristic membrane area 

vector of acceleration of center of mass 

control matrix 

thickness of disk 

gain coefficient of servo valve 

measurement matrix 

moment of Inertia 

membrane stiffness 

constant considering fluid flow losses 

constant considering oil Inertia 

constant considering fluid flow losses 

stiffness coefficient 

stiffness of servo valve 

damping matrix 

shaft diameter 

damping of servo valve 

i dentl ty matr i x 

vector of control forces 

transfer functions (Laplace domain) 

number of degrees of freedom 

external force vector acting on a body 

external force vector acting on a mass element or disk element 
generalized coordinate of servo valve 
gyroscopic matrix 
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G 

h 

I 

dl 

1 

J 

JR 

Jt 

j 

K 

Kpq.fy 

k 

L 

H 

H e 

dM e 

m 

dm 


ffly 

N 

n 

0 

o 

P 


P 

0 


weight 

vector of generalized forces 
tensor of moments of Inertia of a body 
tensor of moments of Inertia of a mass element 
Index; or magnetic bearing current 
performance index 

Jacobian matrix of rotation; 9 u/ 9 q 
Jacobian matrix of translation; 9 r/ 9 q 
index 

stiffness matrix; or gain matrix 

servo valve constants 

gain factor; or number of substructures 

length 

mass matrix 

external torque vector acting on a body 
external torque vector acting on a mass element 
mass 

mass element 

mass Influence of servo valve 
nonconservative matrix 

order of state space representation; or Index 

zero matrix 

Index 

matrix for forces proportional to velocities; or solution of 
algebraic Rlccatl equation; or solution of Ljapunov equation 
pressure 

matrix for forces proportional to displacements; or weighting 
matrix 
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q 

R 

r 

r(z,t) 

s 

T 

T( ju) 

t 

t 

U v 

u 

u(z) 

u(z,t) 

v(z) 

v(z,t> 

X 

X 

X B 

x 0 

y 

y 

z 

a 

0 

6 

X 


vector of generalized coordinates 
outer membrane radius 

Inner membrane radius; or number of actuators 

position vector 

Laplace operator; b + ju 

oil temperature 

transfer function 

torque vector 

time 

control voltage to the servo valve 
control vector 

vector of admissible shape functions 

distributed coordinate (displacement of rotor axis in 

x-di rection) 

vector of admissible shape functions 

distributed coordinate (displacement of rotor axis In 

y-dlrection 

state space vector 

displacement of rotor In x-directlon 
bearing displacement In x-dlrectlon 
characteristic flow velocity 
measurement vector 

displacement of rotor In y-dlrectlon 
axial rotor coordinate 
angle of rotation about x-axls 
angle of rotation about y-axls 

Indicates variation (here only of position); Dirac function 
elgenval ue 
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Sv 

p 

<p 

S<f 

Q 

u 

w 

Uy 

Subscripts 

A 

B 

D 

I 

M 

0 

P 

R 

S 


servo valve damping coefficient 

mass density 

vector of rotation 

small virtual rotlon 

rotor frequency 

vector of angular velocity 

skew-symmetric tensor of angular velocities 

servo valve el genfrequency 

acceleration or actuator 

bearing 

derivative 

Integral ; or Inertial 

magnetic bearing; or measuring location 

output control ler 

proportional 

rotor; or rotating 

sensor location; or state controller 
directions of coordinates 


Superscripts 
e external 

s center of mass 

T transposed 

Mathematical symbols 

a e R n vector a of dimension n x 1 

A e R m * n matrix A of dimension m x n 
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derivative with respect to time; 3/3t 
Integral over entire system 


( ') 
/ 

(S) 

6 


variation with respect to displacements or rotation (no time 
variation) 
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FIGURE 1. - SCHEMATIC OF A MAGNETIC BEARING WITH INTEGRATED SENSORS; 
EHX1, EMX2. EMY1, AND EMY2 ARE ELECTROMAGNETS; GS1 TO 6S4 ARE 
SENSORS. 



FIGURE 2. - SCHEMATIC OF MAGNETIC BEARING INCLUDING SENSOR AND 
ELECTRONIC DEVICES; F x AND F y ARE MAGNETIC BEARING FORCES; u x 
AND U y ARE VOLTAGE INPUTS TO POWER AJFLJFIERS. 
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FIGURE 4 . - ELECTROMAGNETIC ACTUATOR CHARACTERISTICS. 


300 Hz. 


OIL FEED FROM 
SERVO VALVE 



FIGURE 5. - ACTIVE CHAMBER SYSTEM. 



FIGURE 6, - REFERENCE FRAMES AND COORDINATES FOR THE iTH SUB* 
STRUCTURE: 1 INDICATES THE INERTIAL FRAME; R INDICATES THE 
ROTATING REFERENCE FRAME; AND H INDICATES THE BODY FIXED 
FRAME. 


40 


ORIGINAL PAGE IS 
OF POOR QUALITY 





BEARING 2 


ELASTIC 



FIGURE 7, - MECHANICAL MODEL OF ROTOR SYSTEM. 



U) BACKWARD MODES. 


(b) NODES OF NONROTATING 
ROTOR. 


0 = 100 Hz 



E16EN- 

FREOUENCY, 
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Hz 
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82. S 


229. A 
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FIGURE 8. - MODE SHAPES OF THE ROTOR-BEARING SYSTEM FOR ROTOR FREQUENCIES 0 - 0 AND 0 = 100 Hz. 
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(a) RELEVANT NATURAL NODE 
SHAPES, 



(b) HIGH-ORDER NATURAL NODE 
SHAPES. 


FIGURE 10. - GUIDE FOR CHOOSING ACTUATOR AND SENSOR POSITIONS FOR 
MINIMAL SPILLOVER (DOT SHOWS RANGE OF SUCH POSITIONS). 



FIGURE 11. - SIMPLE ROTOR SYSTEM FOR EXPLAINING SPILLOVER 
EFFECTS; AND f Hy ARE ACTUATOR FORCES. 
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FIGURE 13. - MECHANICAL MODEL OF THE EPITAXY CENTRIFUGE. 



150 300 
ROTOR FREQUENCY. 0. Hz 


FIGURE 14. - DAMPING VALUES FOR EIGENVALUES OF EPITAXY 
CENTRIFUGE SYSTEM. 
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FIGURE 15. - ELASTIC ROTOR ACTIVELY IN- 
FLUENCED BY A MAGNETIC BEARING. 


1 FORWARD MODES 

— — — BACKWARD MODES 
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Y~ THIRD 
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Y~ SECOND 


0 100 200 500 kOO 

ROTOR SPEED. 0. H* 

FIGURE 16. - CALCULATED PLOTS OF EIGENFREQUENCIES UP TO FOURTH 
ORDER OF ROTOR SHOWN IN FIGURE 15. CRITICAL SPEEDS OCCUR 
AT IK AND lu. 
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0- 300 Hz 
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- 150 Hz 

0 


EIGENVALUE, 
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(J) NODES OF NONROTAT- <b> FORWARD MODES. <C> BACKWARD MOW 

JNG ROTOR. 

FIGURE 17, - NODE SHAPES OF THE ROTOR SHOWN IN FIGURE 15; M INDICATES ACTUATOR LOCATION. 



(C) BACKWARD NODES. 
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0 100 200 300 *100 500 

ROTOR SPEED, 0, Hz 

FIGURE 19. * DAMPING CURVES FOR OtHPUT FEEDBACK CONTROL. 
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